Single hole dispersion relation for the real CuC>2 plane. 
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Dispersion relation for the Cu02 hole is calculated basing on the generalized t-t'-J model, re- 
cently derived from the three-band one. Numerical ranges for all model parameters, t/J = 2. 4. .2. 7, 
t'/t = 0.. — 0.25, t" It — 0.1. .0.15, and three-site terms 2tN ~ is ~ J/4, have been strongly justified 
previously. Physical reasons for their values are also discussed. Self-consistent Born approximation 
is used for the calculation of the hole dispersion. An excellent agreement between calculated 
and one obtained from the angle-resolved photoemission experiments is found. 
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Recent angle-resolved photoemission experiments 
(ARPES) on insulating copper oxide Sr 2 Cu0 2 Cl 2 [| can 
be considered as a direct test for low-energy models de- 
scribing carriers (holes) in the Cu0 2 plane. Experimen- 
tally observed dispersion relation for a single hole 
has the following characteristic features: (i) bandwidth 
about 2J, (ii) band minimum at the (tv/2, 7r/2) point, (hi) 
isotropic dispersion near the band minimum, and (iv) al- 
most flat dispersion along the line (0,0) — > (ir/2,0) — > 

M). 

The first result agrees with t-J model prediction as well 
as with ones of all possible t-J model generalizations in 
the region of parameters when t > J. Following the phys- 
ical arguments by Kane, Lee and Read Q it seems to be 
rather general that in presence of the strong spin fluctu- 
ations there are no stable quasiparticles at higher ener- 
gies (> 2J). It is of no importance whether the 'bare' 
dispersion of the hole exists (t-t'-J model) or not (t-J 
model) , since the basic arguments are the absence of the 
hole-magnon scattering near the bottom of the band and 
its domination at higher energies. Roughly, it looks like 
some kind of the Cherenkov effect: a massive quasiparti- 
cle cannot create an excitation with linear dispersion up 
to the threshold energy. 

Experimental observation of the band minimum at 
(tt/2, 7t/2) point also agrees with the quasiparticle (spin- 
polaron) dispersion calculated in the framework of t-J- 
like models. It is well established by now that the almost 
degenerate dispersion along the magnetic Brillouin zone 
(MBZ) boundary ((tt, 0) — (0, 7f) line) is an intrinsic prop- 
erty of the pure t-J model, and that it is lifted out by 
any small (compared to t, not J) additional hopping in- 
tegral, for example the next-nearest-neighbor hopping if . 
Moreover, including t'-tcrm in the low-energy model of 
the real Cu0 2 plane is strongly supported by the first 
principle calculations, which show that the direct 0-0 
hopping provide large enough transfer amplitude to the 
next-nearest-neighbor site |3| . 

Thus, it is not much surprising that the masses in the 
directions along and perpendicular to the MBZ boundary 
were found close to each other. To be considered as the 
experimental constraint on the parameters of the t-t'-J 
model at J/t = 0.4 it fixes t' near -0.3t {§. 

Returning to the experimental results, note that theo- 



retical description of the last feature, i.e., flat band along 
(0, 0) — > (tt, 0) line presents a problem. This flat re- 
gion is absent in the t-J model quasiparticle band. Sim- 
ple i'-term adds the 'bare' hole dispersion in the form: 
e k = 4*' cos k x cos k y . Considering this term as the cor- 
rection to the pure t-J model dispersion one can see that 
it does not lift (tt/2, 0) point from its t-J model position, 
and so it cannot provide the flat dispersion. 

Disagreement between experimental band shape and 
theoretical one based on the t-t'-J model returns us to the 
problem of the correct low-energy model of the real Cu0 2 
plane. There were some recent works devoted to this 
problem, which consider Cu0 2 holes in the framework 
of the three-band model in the strong-coupling limit [gj , 
H . These calculations reproduce the experimental band- 
structure much better than the t-t'-J ones, but some of 
the fitting parameters differ from those proposed in the 
cluster analysis of the spectroscopic data and || elec- 
tronic structure works . 

From our point of view the experimental and pure t- 
t'-J model discrepancy is the reason to revise approxi- 
mations were made in obtaining this model for the Cu0 2 
plane, not to deny it. 

In our previous works (^] we developed ideas of the 
three-band model low-energy reduction, firstly proposed 
by Zhang and Rice B. We performed the consistent 
quantitative mapping of the initial model to the single- 
band one using Vannier-ortogonalized basis of the oxy- 
gen states and canonical transformation approach [tLtM . 
It allowed us to obtain the low-energy generalized t-t'-J 
model and to calculate the ranges of its parameters for 
the real Cu0 2 plane pjj ]. 

Our general statement |ll[] is that there are physical 
reasons for including some other terms except t' one, 
namely hopping term to the next-next-nearest neigh- 
bors and the so-called 'three-site' terms, into low-energy 
model. We also should stress that the simple addition of 
the t'-term alone to the t-J model is too naive to give 
the correct description of the subtle details of the hole 
spectrum. 

In this paper we show that the generalized t-t'-J model 
with the set of parameters, which presented in our recent 
work pd| |, reproduces the experimental bandshape at all 
k-points quite well. 



1 



Hamiltonian of the generalized t-t'-J model has the 
form @: 

H = H t _j + H t ,, 

H t -j = t £ Ca<*+ j E s ^> w 

(ij),<x (ij) 
H t > =t' £ L 2 J> + t" 5 i,cAa 

<«i),a (ilj),a/3 

in standard notations of the constrained Fermi operators, 
brackets denote first (()), second and third (()2,3) neigh- 
bor sites, respectively. Three-site terms are written in 
the rotationaly invariant form, Ni, Si are the number of 
fermion and spin operators, respectively. tN and ts ob- 
tained for the usual Hubbard model are 2t^ — —is = 
t 2 /U = J/4, their 'three-band' values are not so simply 
related to the other parameters due to the presence of 
triplet state. 

Ranges for the parameters of the model (Q) arc [ pT| |: 
t/J = 2.4.. 2.7, t! = 0.01.. - 0.25, t" = 0.12. .0.16, t N = 
0.01. .0.07, t s = -0.07.. - 0.16. It is worth to note that 
the t' amplitude is smaller than it follows from the cluster 
calculation Q| and t" is not small compared to t' . As it 
was discussed in Ref. |ll| and in the work by Jefferson et 
al fl^l , the main reason for the difference between CU2O8 
cluster and infinite plane t' hopping parameters is the 
Vannier nature of the latter. It was shown [0, jl2] that 
the Cu-0 and 0-0 hopping amplitudes tend to cancel 
each other for t', and sum up for t" -term. This is the 
cause of not small and weakly waried t" , whereas t" Cu _Q 

(t'6-o)^t'cu-n(.to-o)- 

Previously flll| we have calculated spin-polaron disper- 
sion for the parameters discussed above using the simple 
variational ansatz fll3| ], which is quite good for the pure 
t- J model. It consists of the 'bare' hole and four 'one- 
string' holes, that is, as it is clear by now, not enough for 
the correct treating of the i'-terms. 

In this paper we treat the energy calculation problem 
using the self-consistent Born approximation (SCBA). 
First of all, one should turn to the well known spinless- 
fcrmion Shwinger-boson representation for the Hubbard 
(constrained fermion) operators jL4|. In that case con- 
straint on the fermion degrees of freedom fulfilled exactly 
|jl5|| , and the only approximation is the spin-wave one. 
Hamiltonian of the model ([!]) becomes: 

H ~ Y £ khlh k + Y w q a k a k 

k q 

+ ]T (M k , q 4_ q ft k a+ + H.c.) + (2) 

k,q 

where h^(h), a* (a), are the spinless hole and magnon 
operators, respectively, e k is the 'bare' hole dispersion, 
cj q = 2J(1 — 7 q ) x ^ 2 is the spin-wave energy, M k>q = 



4f(7 k _ q [/ q + 7k Vq), Uq,Vq are the Bogolubov canoni- 
cal transformation parameters. includes the higher- 
order magnon terms. Bare hole dispersion has the form: 

£k = e k + Se k 

el = A{t' + 2t N -t s ){ 1 l-{^) 2 ) (3) 
+ 8(t" + t N - t s /2)( 7k + (7 k ) 2 - 0.5) 

where we used shorthand notations 7k = g(cos(fc x ) + 
cos(fcy)), 7^ = ^(cos(fc l; ) — cos(fc y )). <5ek is the addition 
from zero-point fluctuations: 

5e k = Aax{t' + 2t N + t s ){ll - ( 7k ) 2 ) 

+ 8a 2 (t" + t N + is/2)( 7k + ( 7k f - 0.5) (4) 
+ 4/3i s (4 7 2-l) 

where a\ = 0.138, a 2 = 0.107, [3 — —0.347. One can see 
from Eqs. (^), (^) that in contrast to t! (first terms), t" 
(second terms) lift (n/2,0) point to the higher energy. 

Using SCBA we find the Green function of the hole as 
G(k, u>) = (oj — e k — S(k, uj — ix) q )) _1 with the self-energy 

S(k, W -c q ) = l^A^ q G(k-q,c-c q ) (5) 
q 

It was proved earlier jl6| that the first order correction 
to the hole-magnon vertex is absent and the highest are 
very small. 

Recently Bala, Oles and Zaanen JlTj showed that the 
higher-order vertices (H^) do not change the SCBA re- 
sults and confirmed that one-magnon coupling are accu- 
rate enough to reproduce the realistic properties of the 
t- J-like models. 

Equation (^) was solved numerically by the simple it- 
eration procedure. We found no significant changes of 
the results for 16 x 16 k-points (in MBZ) and 1000 u- 
points, and for 24 x 24 x 3000 points. Also, we checked 
our procedure for the pure t-J model and found very 
close agreement with earlier results fl8|| . Results of our 
generalized t-t'-J model calculations together with the 
experimental points are presented on Fig. |l|. It is im- 
portant to stress that it is not the 'best fit', we simply 
used the average values of parameters from their 'real- 
istic range'. In the main term of the bare dispersion e k 
(^) t j\f and ts terms enter in a combination (2t^ — ts), 
which realistic range is [|J.|J], so we simply take its 
Hubbard value (2t N - t s ) = J/2 = 0.2t. We used 
t/J = 2.5, t N = J/8,< s = - J/4, and J = 0.14 eV (from 
Ref. H). For t! and t" we used t! = -0.2t, t" = 0.15i. 
Note, that the bare dispersion (|3|) consist of two terms, 
which can be considered as the t' e ^ (= t' + 2t^ — ts) 
and t'/ff (= t" + t N — is/2), and since t' and (2t^ — t s ) 
have the opposite signs, t' e ^ becomes very small at all 
realistic values of parameters. 

Figure |l| shows an excellent agreement with experi- 
ment along (it, 0) — (0, n) as well as along (0, 0) — (n, 0) 
lines. Notice, that another feature of the ARPES can be 
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explained in the spin-polaron approach. Lower intensity 
of the photoemission peaks at the top of the hole band 
(bottom in the electron language) easily connected to the 
lower quasiparticle residue at (0,0) — > (71", 0) points. The 
last is due to importance of multi-magnon scattering pro- 
cesses for the 'cutting' of the wide initial ('bare') band 

We also found an important feature of the energy spec- 
trum of the proposed model: if the values of t", ts, *jv 
are not small (average and larger) , the shape of the band- 
width is fully insensitive to the t' changes. Changes of 
*' only shift the energy of ground state and change the 
quasiparticle residue at the top of the band. Opposite to 
it, t" strongly varies the (ir/2,0) position. These results 
are shown in Fig. |^. They are easily understood remem- 
bering that the higher energy states are unstable. When 
t" and three-site terms are not very small, they already 
form the band up to the characteristic energy 2J, and 
further changes of t' (even in a broad region) touch only 
the higher states, which are unstable. It is interesting 
that the further increasing of t" (> 0.15*) also do not 
change the shape of the band. 

Summarizing, we showed that the generalized version 
of the t-t'-J model accurately derived from the three- 
band model describe very well the experimental results 
of ARPES on Sr 2 Cu0 2 Cl2 system §. Parameters of the 
model, at which good agreement achieved, are from real- 
istic regions 11 and so they are strongly justified. Hop- 
ping integral to the next-next-nearest neighbor site (*") 
as well as three-site hopping terms (ts, *jv) is found to be 
the key parameters for the description of the flat region 
along (0, 0) — (it, 0) line. It is argued that the isotropy of 
the spectra around the minimum (tt/2,tt/2) easily arise 
at any (not very small) t', t", ts, *jv parameters of the 
definite sign. In addition we found that the shape of the 
spectrum is insensitive to varying of t! if the other param- 
eters are not small. Thus, the model has some rigidity to 
parameter changes. Smaller intensity of the photoemis- 
sion peaks at the top of the band can be directly related 
to the small quasiparticle residues at these points. 
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Figure captions 



FIG. 1. Dispersion curve of a hole in the gen- 

eralized t-t'-J model (|l|), (Q) along the main directions 
(0,0) -» (tt.tt), (tt,0) -» (0,0), and (tt,0) -> (0, tt) 
(solid curve). Model parameters that provide this i?k are: 
t/J = 2.5, t' = -0.2*, t" = 0.15*, t s = -2t N = - J/4, J = 0.14 
eV. Experimental results from Ref. [jjj are also shown (open 
circles) . 



FIG. 2. Dispersion relation of a hole in the gen- 
eralized t-t'-J model for the different sets of parameters: 
t' = -0.2t,t" = 0.15* (solid curve), t' = -0.5*,*" = 0.15* 
(dashed curve), *' = —0.2*,*" = (dotted-dashed curve). 
Other parameters are: t/J = 2.5, *s = — 2*jv = — J/4. 
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